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SUMMARY 

In this paper we compare the results of Donato's exponentially varying ground model, 
Attenborough's exponentially varying ground model and the rigid backed thin layer model. We 
show that these models produce similar results for slow variations. For rapid variations the results 
are quite different but the basic theory used is only correct for the thin layer model. These results 
suggest that the exponentially varying models are not necessary for fitting ground impedance data. 

INTRODUCTION 

Donato proposed an exponentially varying ground model to be used for the interpretation of 
ground impedance data.* Attenborough has demonstrated that the exponential variation chosen by 
Donato results in model grounds with increasing porosity with depth and has derived a ground 
model which has a decreasing porosity with depth.^ 

In this paper we examine the behavior of both these models in the limit of large and small 
variation and compare the results to the rigid backed layer model. 3 To facilitate this we have 
reduced the solutions to their simplest forms and have employed Attenborough's low frequency/high 
flow resistivity results for numerical comparison. 

I. GROUND MODELS 


A . Rigid Backed Layer 

A layer of porous material of thickness d overlying an acoustically rigid surface has a surface 
impedance of the form: 


Z(0) = i Z c cot (kd) (1) 

where Z c is the impedance of a seminfinite half space of the porous material and k is the complex 
wave number in the porous material. 
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B . Donato's Exponential Model 


Donato has derived a impedance model for a material whose porosity times wave number 
decreases exponentially with depth. Attenborough has demonstrated that for natural grounds this 
implies that the porosity increases exponentially with depth and the wave number decreases 
exponentially with depth. This will not commonly occur in natural ground surfaces but may be a 
useful model in some cases. With the notation above Donato's formula becomes 


Z(0) = iZ c 


Jo(2k/a) 
Ji(2k/a) ’ 


( 2 ) 


a is the exponential varation of the square of the complex wave number 

k(z) 2 = k(0) 2 e- az . 


(3) 


C. Attenborough's Exponential Model 

Attenborough's solution for a porous material whose porosity decreases exponentially with 
depth and wave number increases exponentially with depth is given by 


Z(0) = iZ c 


Hj > 2) (2k/a) 
Hj 2) (2k/a) ’ 


(4) 


where 


k(z) 2 - k(0) 2 e az . (5) 

II. BEHAVIOR OF THE IMPEDANCE AND WAVE NUMBER 

It will be useful in the interpretation of these models to have a specific formulae for the wave 
number and impedance of a homogeneous porous material. For this paper we will use 
Attenborough's low frequency approximation: 

Z C = J&- =.218(^-Y /2 (1 + i). (6) 

yQco Vf/ 

o e is the effective flow resistivity of the material, y is the ratio of specific heats and c is the speed of 
sound in air. 2 
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m. BEHAVIOR OF THE GROUND MODELS IN THE LIMIT OF LARGE AND 

SMALL ARGUMENTS 


A . Rigid Backed Layer 

i) Limit as d — » 0. 

For a thin layer d — » 0 and Eq. (1) becomes 

(7) 


( 8 ) 

where is co/c. Note that the imaginary term approaches infinity as k Q d goes to zero, while the 
real part depends only on the layer thickness and the surface flow resistance. This form is 
displayed by Attenborough. 2 

ii) Limit as d — > oo 

As d — » °° the model should recover the result for the homogeneous half space. The 
cotangent can be expanded in terms of the exponents of the real and imaginary parts of kd. 

lim cot(kd) = lim = e:*' d e+M / eijElLejM.,- s-*» d e*** (9) 

d— d— 2 2i 


Z(0)= lim iZcCot(kd) = i^ -i^ 
d— >0 kd 3 

If we use Eq. (6) to relate Z c and k for low frequency we find 

Z (0)- 4,I (- 218 ,) 27Qdff %i-J— 

30 -ySlkod 


where 


k = kj + i k 2 . 
k 2 must be positive so that 


Z(0) = iZ c (-i) = Z c , 


GO) 


and the original condition is recovered. 
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B . Donato's Exponential Model 

i) Limit as a becomes small 


As a becomes small the medium approaches a homogeneous media. If we take the limit of 
Eq. (2) for small a and large 2k/a we find 


Z(0) - i Z c 


.J 2k<x cos /2k . K ^ 

V 2k V a 4 / 

J 27ta rn o /2k. 7t _7t \ 
V 2k ° \a 2 4) 


= i Z c cot 


/ 2k 7L \ 
Va 4 / 


(ID 


This is like the impedance of a thin layer of thickness 2/a with an additional -it/4 phase 
change. The next correction term is of order a/2k. A pressure release backed thin 
layer would have a phase change of -ic/2. As a — > 0, the cotangent term will approach -i as in 
Section A-ii) and Z(0) = Z c as expected. 


ii) Limit as a — » 

In the limit as a — > the argument becomes small and the ascending series may be used to 
evaluate the Bessel functions. 


zm “ iz <>i:- i -te • <12) 

a 

The behavior of this solution is very similar to Eq. 7. We have a rapidly increasing imaginary 
part and a constant real part as the frequency decreases for fixed d and o e . The imaginary parts are 
identical if the rigid backed layer has a thickness 1/a, while the real parts are equal if the rigid 
backed layer thickness is given by 1.5/a. 

C. Attenborough's Exponential Model 
i) Limit as a — » 0 

The asymptotic expansions can be employed for the Hankel functions giving 

. / 2 a q — i( 2 kJa - 71/4) 

Z(0) - i z c Yrc& = i z c e- i7t/2 = z c . (13) 

./2q e -i(2k/a - nU) 

» it 2k 62 6 
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The Attenborough model recovers the homogeneous half space surface impedance as a — » 0. 
ii) Limit as a — » CO 


The small argument formulae for the Hankel functions are inserted in Eq. (4) to give 


Z(0) = Z c [(|-ie)^-i^ln(|)‘, 


(14) 


where e = .5772. 


This result is not easily interpreted in terms of a layered model. The behavior of this solution 
is best illustrated by use of Eq. (6) to yield 


Z(0) = 5.923 In 


+ 3.419 + i 13.955 


(15) 


As a — » <x> the impedance of the Attenborough model has a large negative real part tending to - «> 
and a constant imaginary part This puzzling result indicates that the surface is not absorbing 
energy and has a reflection coefficient greater than one! In a gross sense the behavior is physical. 
The reflection coefficient approaches one as the impedance becomes infinite. The only problem is 
that the surface cannot be generating acoustic energy. 


iii) Limit for 2k/a > 1, a not infinite 

A third limit is developed by Attenborough as useful for computation and comparison with 
data. This form is developed for a small enough that the leading term in asymptotic series for the 
Hankel functions may be used. For 2k/a > 1 


H< 2) (2k/a) , Q+jf) . . r, , B] 

H ( , 2) (2k/a) (l ' 8 ^ 2k) 4k '* 


(16) 


and 


Z < 0) * Z c{ 1+ tk} 


Using Eq. (6) to relate k and Z c gives us 


Z(0) = Z c + 


ic Ta 

4 Y®La ‘ 


(17) 


(18) 
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Defining a e = a/Q and inserting numerical values from Eq. (6) gives us Attenborough's form: 

Z(0) = .218 (y) 1/2 + i .21 8 (y - ) 1/2 + 9 - 74 (y)] 0 9 ) 

The next terms in the asymptotic series are on the order of 7% of the last term in Eq. (19) when the 
argument of the Hankel function is one. 

Note that we can recover Eq. (13) by letting a approach zero. Also note that the second term 
in Eq. (17) is very similar to the form for the imaginary part of the impedance of a thin rigid backed 
layer. Compare 

iZ c -SL and %■ . (20) 

c 4k kd 

The second term in Eq. (17) is the imaginary part of the impedance of a thin layer of effective 
thickness d e = 4/a. The imaginary parts dominate the impedance for large a. 

IV. NUMERICAL RESULTS 

To calculate numerical values for the three impedance models we set 


kd = 2k/a = x(l + i). 

Then, using Eq. (6), we solve for f and Z c in terms of x: 

p __ axe l 2 1 
[47cyQ(.218)J ’ 

and 


( 21 ) 


( 22 ) 


47ryD(.218) 2 o e 


(1 + i) 


We use the following typical values of y, Q, o e , and a based on our experience and that of 
Attenborough: 
y = 1.4 
Q = 0.4 

a e = 120,000 MKS rayls 
a = 40. m' 1 ; d = 5 cm. 

Then, we calculate impedances using Eqs. (1), (2) and (4) for x = 0 to 5. The results are plotted in 
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Fig. 1 (rigid backed layer), Fig. 2 (Donato's formula), and Fig. 3 (Attenborough's solution). The 
imaginary parts of the impedance are multiplied by -1 so the plots of the real part are usually on the 
positive side of the vertical axis and the imaginary parts are on the negative side. The plots are 
nearly identical for values of x greater than one. For the variables above, x = 1.0 corresponds to 
654 Hz. 

Figure 4 displays the normal reflection coefficient calculated from Eqs. (1,2, and 4). The 
behavior is similar for all the models. Better agreement can be achieved between any two models 
by the choice of the equivalent depth of the exponential variation. 

V. DISCUSSIONS AND CONCLUSIONS 

The surface impedance predicted by each of the three models above approaches the 
homogeneous half-space impedance as the variation of wave number becomes small or the layer 
depth becomes large in the rigid backed model. 

As the exponential variations become larger the impedance formula can be approximated as a 
constant or slowly varying real and imaginary part plus an imaginary term which is proportional to 
cx/0) or 1/oxl. 

For very rapid variations, the expansion of Attenborough’s solution results in a non-physical 
solution (Eq. 13). 

The basic assumption in the derivation of Eq. (6) and it's more exact analogues, is that the 
gradients of the variables with respect to the propagation direction are much smaller than gradients 
of the variables normal to the direction of propagation.^ The result that the reflection coefficient is 
greater than one for small variable x is probably due to the error in Eq. (6) rather than any physical 
error in the theory leading to Eq. (4). 

By the same reasoning, Donato's formula should be inaccurate for small values of the variable 
x. There is no physical problem with the thin rigid backed layer since the porous layer is 
homogeneous and Eq. (6) should hold. For the variables we have chosen, there appears to be little 
practical reason to employ the exponential models to fit ground data, while there appears to be a 
significant theoretical reason for not using the exponential models in the region where they vary 
significantly from the rigid backed layer. 

At very low frequencies, the impedance translation theorem can be employed to calculate the 
impedance of an impedance backed layer. This model has sufficient flexibility to fit most data 
without the theoretical difficulties of the Donato or Attenborough models. 
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Figure 1. The real and imaginary part of the grbund impedance versus the parameter x for the 
thin rigid backed layer. The imaginary part is multiplied by negative one for 
display purposes. 
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Figure 2. The real and imaginary parts of the ground impedance versus the parameter x for 
Donato* s exponentially varying model. 
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Figure 3. The real and imaginary parts of the ground impedance versus the parameter x for 
Attenborough’s exponentially varying model. 



Figure 4. The absolute value of the reflection factor for normal reflection for the three 
models versus the parameter x — thin layer, Donato, Attenborough.) 




